	*===========================*
	* Table  - suspensions
	*===========================*
	preserve
		
			* look at stsusp and ltsusp and exp
			tab subject, gen(__tmp_subj)
			
		*quiet{
			label variable log_ltsusp "Long-term Suspensions (Log)"
			label variable log_stsusp "Short-term Suspensions (Log)"
			label variable black "Black"
			label variable hisp "Hispanic"
			label variable white "White"
			
			
			label var exp_capped "Experience"
			label var exp_capped_squared "Experience Squared"
		
		
			global subjectvar __tmp_subj1 __tmp_subj2
			
			drop subject
			duplicates drop stu_id sch_id id, force
			eststo clear
			estimates clear
			
			* generate insample
			reghdfe stsusp $experiencevar $nothreatvars $subjectvar $yeardosevar, absorb(grade stu_id sch_id id year) vce(cluster id#year)
			gen _insampletmp = e(sample)
			keep if _insampletmp
			local j = 1
			estimates clear
			eststo clear
			
			reghdfe stsusp $nothreatvars $subjectvar $yeardosevar, absorb(grade stu_id sch_id year) vce(cluster id#year)
			estimates store model_`j'
			distinct id if e(sample)
			estadd scalar _ndistinct = r(ndistinct)
			quietly: sum stsusp if e(sample), meanonly
			estadd scalar _meanoutcome = r(mean)
			local j = `j' + 1
			
			reghdfe stsusp $nothreatvars $subjectvar $yeardosevar, absorb(grade stu_id sch_id id year) vce(cluster id#year)
			estimates store model_`j'
			distinct id if e(sample)
			estadd scalar _ndistinct = r(ndistinct)
			quietly: sum stsusp if e(sample), meanonly
			estadd scalar _meanoutcome = r(mean)
			local j = `j' + 1
					
			reghdfe stsusp $nothreatvars $experiencevar $subjectvar $yeardosevar, absorb(grade stu_id sch_id id year) vce(cluster id#year)
			estimates store model_`j'
			distinct id if e(sample)
			estadd scalar _ndistinct = r(ndistinct)
			quietly: sum stsusp if e(sample), meanonly
			estadd scalar _meanoutcome = r(mean)
			local j = `j' + 1
			
								
			* prepare estimates for estout
			estfe model_*, labels(sch_id "School FE" year "Year FE" grade "Grade FE" stu_id "Student FE" id "Teacher FE")
			return list
		
		
		esttab model_* using "$output_dir\Tables\_Table_18.tex", replace ///
			indicate(`r(indicate_fe)' "Experience = *xp*", labels("X" " ")) ///
			keep(nothreat_cycle3)  noomitted nobase se label ///
			nonotes $stars substitute(/ "") starlevels(\$^+\$ 0.10 * 0.05 ** 0.01 *** 0.001) ///
			scalars("_meanoutcome Average Suspensions" "_ndistinct Teachers") noobs
		* return stored estimates to their previous state
		estfe model_*, restore
		
		
		
		
	restore
